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The problem of both electron and phonon scattering by wedge disclination dipoles (WDD) is 
studied in the framework of the deformation potential approach. The exact analytical results for 
the mean free path are obtained within the Born approximation. The WDD-induced contribution 
to the residual resistivity in nanocrystalline metals is estimated. Using the WDD-based model of 
a grain boundary, the thermal conductivity, «, of polycrystals and amorphous dielectrics is stud- 
ied. It is shown that the low-temperature crossover of k experimentally observed in LiF, NaCl, 
and sapphire can be explained by the grain-boundary phonon scattering. A combination of two 
scattering processes, the phonon scattering due to biaxial WDD and the Rayleigh-type scattering, 
is suggested to be of importance in amorphous dielectrics. Our results are in a good agreement with 
the experimentally observed k in a-SiC>2, a-Ge02 a-Se, and polystyrene over a wide temperature 
range. 
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I. INTRODUCTION 

There exist many varieties of extended defects in crys- 
tals, topological in their origin. The most known exam- 
ples are dislocations, disclinations, twins, grain bound- 
aries, stacking faults etc. These defects play a significant 
role in description of various phenomena in real crys- 
tals as well as in non-crystalline materials. In particular, 
there is reason to believe that linear defects like dislo- 
cations and disclinations are the principal imperfections 
of liquid crystals!, some amorphous solidsBEl, and poly- 
mersD. 

The contribution to the transport characteristics 
due in. dislocations is now well understood (see, e.g., 
RcfsB'El). Some aspects of the qualitative behavior of 
the disclination-iaduced electron scattering have been re- 
cently presented!! In particular, it was found that both 
dislocations and disclinations can be effective scatter- 
ing centers for conducting electrons, especially at low 
temperatures when other scattering mechanisms are sup- 
pressed. Thus, along with point impurities, these defects 
give a contribution to the residual resistivity. In real crys- 
tals, however, the isolated disclinations are rather exotic 
objects. Instead, for topological reasons, the dipole con- 
figurations are more favourable. In addition, the creation 
energy for a single disclination considerably exceeds that 
for a disclination dipole. 

There is another reason to call attention to dipoles 
of disclinations. As a matter of fact, the wedge discli- 
nation dipoles (WDD) simulate finite dislocation walls. 
In turn, dislocation walls describe the low-angle grain 
boundaries!'!! Thus, one can expect that the results ob- 
tained for disclination dipoles will be useful in description 
of the grain-boundary scattering problem. This allows 



us to extend the possible applications to a study of the 
transport properties of polycrystals where grain bound- 
aries are of importance. In particular, there is experimen- 
tal evidence that the grain boundaries give contribution 
to the resistivity in metals (see, e.g., review in Ref.EEl) 
which in part depends on the size of the grain boundary. 
Notice that though the problem of the grain-boundary- 
induced scattering has been formulated many years ago 
the proper solution is still absent. 

A theory of the phonon scattering by grain boundaries 
has [been developed within the Born approximation in 
Ref.O. A grain boundary was considered as a wall of edge 
dislocations with a rather strong assumption that the dis- 
location wall is infinitely long. Nevertheless, the main 
properties predicted by this model have been confirmed 
by experiments£3~li-l In particular, the phonon mean free 
path was found to be a constant over a wide temperature 
range. As a result, the low-temperature thermal conduc- 
tivity in polycrystals varies as T 3 in agreement with the 
experimental data. However, the model of an infinitely 
long dislocation wall failed to describe a remarkable in- 
crease in thermal conductivity observed below-scpae char- 
acteristic temperature T (T ~ 0.1K for LiPlMil). 

Recently!!!, the grain boundary phonon scattering 
problem has been investigated within the more realis- 
tic model which takes into account the finiteness of the 
boundary. The basis for this model was the known anal- 
ogy between-diaclination dipoles and finite walls of edge 
dislocationstSIIZl. It was found that the proper consid- 
eration of the finiteness of the boundary results in the 
low-temperature crossover of the, thermal conductivity in 
agreement with experimentsli^O. 

A more intriguing application is an attempt to de- 
scribe the physics of amorphous dielectrics by considering 
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WDD as the basic structural elementstS. This point of 
view agsces with the cluster picture proposed earlier for 
glassestaE3. In accordance with this picture there exist 
many small crystalline grains (microclusters with average 
diameters of order 20-30A) in the amorphous state. On 
the other hand, the concept of elastic dipoles for orienta- 
tional glasses was introduced and successfully explored in 
RefsElo. The authors considered elastic dipoles as the 
additional to two-level systems (TLS) structural elements 
of glasses. This assumption allows them to describe both 
the specific heat and the thermal conductivity of some 
glasses in the plateau region and above (up to 100K). 
The physical nature of these dipoles, however, has not 
yet been clarified. As is knownH, a possible way to un- 
derstand the microstructure as well as the nature of prin- 
cipal imperfections is to study the transport properties 
of a material. It will be shown below that the concept of 
the elastic dipole can be successfully realized via dipoles 
of wedge disclinations without considering TLS. Namely, 
the WDD-induced phonon scattering is found to provide 
the correct description of the low-temperature thermal 
conductivity of amorphous dielectrics. Notice also that in 
accordance with the geometrical consideration disclina- 
tionsj-are expected toralways be present in the amorphous 
stateEa (see also Reffl). 

The main goal of the paper is twofold. First, we out- 
line the general formalism of the WDD-induced scatter- 
ing for both electrons and phonons. The most important 
details needed for a better understanding of calculations 
are given. Second, we apply the results obtained for de- 
scription of two important problems. The first one is the 
experimentally observed deviation of the thermal con- 
ductivity from a T 3 dependence below 0.1 K in LiF and 
NaCl. The second problem is the thermal conductiv- 
ity of amorphous dielectrics. It was shown in Rcf.li3 that 
the experimentally observed behavior of thermal conduc- 
tivity of a-SiC>2 over a wide temperature range can be 
explained by a combination of two scattering processes. 
The first one comes from the phonon scattering due to 
biaxial WDD while the second one is the Rayleigh type 
scattering. In the present paper we extend our investiga- 
tion to other glasses and discuss some unsolved questions. 

The paper is organized as follows. The general for- 
malism of the WDD-induced scattering is presented in 
Sec. II. We consider all possible types of WDD and cal- 
culate the corresponding phonon mean free paths. The 
principal distinction between the scattering properties of 
the uniaxial and biaxial dipoles is shown. The developed 
approach is applied to estimation of the WDD-induced 
contribution to the residual resistivity of granular metals 
in subsection A. The phonon scattering due to WDD is 
studied and the contribution to the thermal conductivity 
is calculated in subsection B. The obtained results are 
compared with the experimental data for LiF. In Sec. Ill 
we apply the WDD-induced mechanism of phonon scat- 
tering to the problem of thermal conductivity in amor- 
phous dielectrics. The results are compared with the 
experimentally observed k in a-Si02, a-GeC>2 a-Se, and 



polystyrene (PS). Sec. IV is devoted to the detailed dis- 
cussion of the results obtained, specifically with relation 
to the proposed WDD-based model for dielectric glasses. 



II. THEORY OF THE WDD-INDUCED 
SCATTERING 

In this section we study the contribution to the ef- 
fective cross-section which comes from the potential as- 
sociated with a static deformation of a lattice caused 
by straight WDD. Two reasonable approximations are 
in common usage in studies of such problems (see, e.g., 
Refs.Bfl). First, we suppose that the scattering processes 
are elastic and, second, the Born approximation is valid. 
Besides, we will consider here the simplest case when the 
only elastic deformations are dilatations. 

In this case, an effective perturbation energy due to 
the strain field caused by a single WDD takes the form 

U(r)=GSpE iji (1) 

where Eij is the strain tensor and G is an interaction 
constant. 

Let the disclination lines be oriented along the z-axis, 
the position of the positive disclination be (0,-L) while of 
the negative one (0,L) (see Fig.l). 
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FIG. 1. Schematic picture of a wedge disclination dipole 
with shifted by distances h and h axes of rotation. The 
disclination lines are oriented along the z-axis for both de- 
fects 

Notice that in the general case the axes of rotations 
(with Jli = £le z and = — ^e z ) are shifted relative 
to their positions by arbitrary distances l\ and 1%, re- 
spectively. For lx = l<z = one gets the biaxial WDD 
with nonskew axes of rotation. It was showrO that this 
dipole can be simulated by a finite wall of edge dislo- 
cations with parallel Burgers vectors. In particular, the 
far strain fields caused by biaxial WDD agree with those 
from a finite dislocation wallEZl. 

For l\ — I2 = 2L and li = —I2 one gets the uniaxial 
and the symmetrical uniaxial WDD, respectively. Notice 
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that uniaxial WDD can be simulated by a finite wall of 
edge dislocations complemented by two additional edge 
dislocations at both ends of the wall. The sign of these 
dislocations is opposite to that of the dislocations in the 
wall and absolute values of Burgers vectors are equal to 
b = 2itan(n/2) (b = b y at chosen geometry). As a 
result, the uniaxial WDD becomes a strongly screened 
system as opposed to the biaxial WDD (see Fig. 2). In 
the general case, h ^= h i= 0, one gets the biaxial WDD 
with shifted axes of rotation. 
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FIG. 2. Contour plot of the perturbation energy (2) 
U(x,y)/B in arbitrary units for biaxial WDD (h = h — 0, 
upper picture) and for uniaxial WDD (h — h — 2L, lower 
picture). 

Substituting the explicit form of Eij into Eq.(l) (see 



Apendix A) we find for the perturbation energy U(r) 

x + L 



U(x,y)=B 



(x-L) 
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(x + L) 2 + y 2 



+ k 



x — L 



(x - L) 2 + y 2 



(2) 



where B = Gv(l - 2ct)/(1 - a), v = D,/2tt is the Frank 
index, and a is the Poisson constant. Notice that all 
possible types of WDD are included in Eq.(2). 

As is seen from Eqs.(Al) and (A2) in Appendix, all 
strains caused by WDD are located in the xy-planc. In 
this case, the only components of the wavevector that 
are normal to the defect lines, q = q^, are involved in 
the scattering process. For the sake of simplicity, let us 
assume that incident carriers are normal to disclination 
lines. 

The problem reduces to the two-dimensional case with 
the mean free path given by 



A" 1 =m (1 - cos 9)U(9)d9. 
Jo 



(3) 



Here 5ft(#) is an effective differential scattering radius, 
and rii is the areal density of WDD. Notice that when 
axes of WDD are oriented randomly, one has to perform 
an additional averaging over all possible angles of inci- 
dents. As was shown for dislocationsD, however, such av- 
eraging leads merely to a modification of the numerical 
factor in K(6>). 
^Vithin the Born approximation, K(#) is determined to 
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\<q\U(v)\q > 



(4) 



where all vectors are two-dimensional, 5* is a projected 
area, v is a carrier velocity, the bar denotes an aver- 
aging procedure over a which defines an angle between 
p = q — q and the ir-axis. In other words, it means 
the averaging over randomly oriented dipoles in the xy 
plane. Evidently, the problem reduces to the estimation 
of the matrix element in Eq. (4) with the potential from 
Eq.(2). For this purpose, it is convenient to use the polar 
coordinates (r, <p) 

U{p,a) = <q\U(r)\q > 

= i J d 2 r exp[zprcos(</> — a)]U(r 7 (j)). (5) 

For elastic scattering, the matrix element in Eq.(5) de- 
pends only on |q| = |q | and the scattering angle 9. Thus, 
p=|p| = |q-q|=2gsin(0/2). 

After integration in Eq.(5) and following averaging of 
\U(p,a)\ 2 with respect to a one obtains (see details in 
Appendix B) 
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where A; = h — h, J n {t) are the Bessel functions. Upon 
integrating Eq.(3) with respect to one finally obtains 



, B 2 L 2 TT 2 ni f 2 /l 2/ ,\ 4 , 
= ) U + J o(^L)y — J Q (2 q L)J 1 (2qL) 



^ (j 2 (2gi) + J 2 (2 g L)) j, (7) 



where z = Ai/L. It should be emphasized that Eq.(7) is 
the exact result which allows us to describe all types of 
WDD. Notice, that the behavior of A in Eq.(7) is actually 
governed by the only parameter 2L. 

Let us consider two important limiting cases. For bi- 
axial dipoles with A; — (z = 0) one obtains 



A^ = 8g2 S 7r2 {^o 2 (2^) + J 1 2 (2 g X) 
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M2qL)J 1 (2qL) 



(8) 



For uniaxial dipoles A; = 2L (z = 2), and the mean free 
path A is 



, , 2B 2 L 2 n l -K 2 „ , 
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1 + J£(2qL)- J({2qL) 



- — J (2qL)J 1 (2qL) [ >. (9) 



In what follows, we apply the developed formalism to 
the problem of the WDD-induced electron and phonon 
scattering. 



A. Electron scattering: WDD-induced residual 
resistivity in metals 

As is known!, the residual resistivity of metals may 
be caused by electron scattering due to linear defects 
like dislocations, stacking faults and grain boundaries. It 
was mentioned in the introduction that the WDD-based 
model is a good candidate for modelling the grain bound- 
aries. Thus, the previous analysis allows us to study 
the contribution to the residual resistivity due to grain 
boundaries. 

Let us use the well-known Drude formula for the resis- 
tivity in the static regime 



P = ( — )(-), 
ne A t 



(10) 



where r is the relaxation time, m and e are the mass and 
(5) the charge of the conducting electron, and n is the elec- 
tron density. For point impurities and linear defects like 
dislocations and disclinations, the angle brackets denote 
the configurational average. In our case, this is averaging 
over a in Eq.(4). Thus, one can write the final result 



^(^f)A- 1 



(11) 



with A from Eq.(7), where q = qF, v = vp, and G = Gd- 
Obviously, index F denotes the Fermi values, and Gd — 
(2/3)Ep is the deformation potential constants, where 
Ep is the Fermi energy. For simplicity, we restrict our 
consideration to metals of zinkblende or wurtzite struc- 
tures. In metals typically qp w 0.6 ~ 1.2A -1 . A dipole 
separation is chosen to be 2L w 10 2 ~ 10 4 A which is 
of order of the grain size in polycrystals. In this case, 
2qpL >> 1 and, consequently, the Bessel functions in 
Eq.(7) can be approximated by their asymptotic values. 
The results are the following 
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(13) 



where B e = B(G = G d ). 

The main difference in the behavior of pu and p U n% 
comes from their L-dependence. As is seen, at fixed m 
Pbi ~ L while p un i ~ L 2 . As the result, pu is found to be 
larger than p un i- This agrees with the above-mentioned 
properties of these dipoles. In the case of biaxial WDD, 
the main contribution comes from the low-angle scat- 
tering processes in view of the long-range character of 
the perturbation energy in Eq.(2). On the contrary, 
for uniaxial WDD the large-angle scattering dominates 
since they are strongly screened systems (see Fig. 2). In 
particular, at 2L - 10 3 10 4 A, n, ~ 10 9 cm~ 2 , and 
v = 0.01 one obtains pu ~ 10~ 12 -j- 10~ n il cm while 
Puni ~ 10~ 9 -r 10 -7 cm. Notice that the value of p un i 
agrees witL-thajt in the case of edge dislocations of a sim- 
ilar densityE3'E3. 

It should be noted that the experiments show an in- 
crease of the residual resistivity of nanocrystalline metals 
with L decreasing!!]. This can be explained within the 
above-proposed model due to a direct L-dependence of 
m. Indeed, it is reasonable to assume that rii ~ L~ p 
with 1 < p < 2 (usually p = 2). Thus, in accordance 
with Eq.(12) pu increases with L decreasing. It is in- 
teresting to note that in this case p un i decreases with 
L (except p = 2), as is seen from Eq.(13). Let us reit- 
erate that only the biaxial WDD with nonskew axes of 
rotations (BWDD) simulates the grain boundary. 
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B. Phonon scattering: low-temperature heat 
transport in polycrystals 



It is clear that phonon scattering due to WDD will alsc 
affect the low-temperature thermal conductivity, k. We 
start from the well known kinetic formula 

K= - / C(uj,T)v s l{cj,T)dcj, (14) 
Jo 



where C(uj,T)duj is the specific heat contributed bj 
acoustic phonons within the frequency interval duj, v t 
is an average phonon velocity, and 1(uj,T) is the phonor 
mean free path. It is suggested that C(lu,T) has stan- 
dard form with quadratic in ui density of states and the 
Debye cutoff too- 

The effective perturbation energy caused by the strair 
field of WDD is determined asBB U(r) = hui-/SpE AB . 
where Tito is the phonon energy with wavevector q. 
uj = qv s , v s is the sound velocity (for simplicity it it 
assumed that three acoustic branches are equivalent) 
and 7 is the Griineisen constant. As previously, we sup- 
pose that incident phonons are normal to the disclina- 
tion lines, so that we deal with the two-dimensional scat- 
tering problem. The principal difference from the case 
of electron scattering is the explicit g-dependence of the 
perturbation energy (see details in Reffl). Namely, the 
strain tensor due to WDD remains the same while the 
coefficient B in Eq.(2) should be replaced by B p h — 
hqv s "fv(\ — 2a)/ (1 — a). Eqs.(4),(5) preserve their form 
in this case as well. 



For BWDD the phonon mean free path is found to b< 




/- 1 = 2D 2 (2vL) 2 n iq [j 2 {2qL) + jf{2qL) 

-^Z J o(2<lL)Ji(2qL)), (15) 

where D = tt^(1 — 2cr)/(l — a). For uniaxial WDD one 
obtains 

C!i = ^-(2uL) a mq(l + J 2 (2qL) - jf{2qL) - 

■^J (2qL)J 1 (2qL)y (16) 

Fig. 3 shows 1(uj) for three types of WDD. 
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FIG. 3. Phonon mean free path l(uj) at 2qoL = 6 x 10 3 
for A; = (solid line), A; = 2L (dotted line), and 
A; = 0.5L (dashed line). The parameter set used is: 
L = 1.35 x W' 5 cm, v = 0.023, K = 2.6, m = 1.8 x 10 7 cm -2 , 
and v s = 4.8 x 10 5 cm/ sec. 

We have used a typical for polycrystals size of the 
grain boundary 2L = 2700A. As is seen, three curves 
behave differently. At low frequencies the scattering 
by uniaxial WDD resembles that by a point impurity. 
Namely, it strongly depends on oj, l un i ~ a; -5 , thus once 
again confirming a view of uniaxial WDD as a strongly 
screened system. At high frequencies uniaxial dipoles 
scatter phonons like dislocations with l un i ~ lo^ 1 . It is 
interesting that the same w _1 -dependence appears for ar- 
bitrary biaxial dipoles both at low and high frequencies. 
What is more important, there is the only type of biaxial 
dipoles, BWDD, which shows the unique behavior with. 
Im — » const as w increases (see Fig. 3). It was foundE£l 
that the change in behavior of lu occurs at 2qL ~ 1 or, 
equivalently, at lu* ~ v s /2L. It should be emphasized 
that this intriguing result provides the basis for the fol- 
lowing important speculations. Notice that some visible 
irregularities in Fig. 3 came from rapid oscillations of the 
Bessel functions near the characteristic frequency to*. 
The BWDD-based model was successfully usedibr de- 
scription of the phonon transport in polycrystalstJ. For 
this purpose, Eq.(14) should be integrated with the 
phonon mean free path from Eq.(15). The result is shown 
in Fig. 4 together with the experimental data for LiFt3. 
As is seen, the predictions of the theory are in good agree- 
ment with the experimental results. It is interesting to 
note that the behavior of k(T) in our model is governed 
mainly by D, 2uL and in Eq.(15), that is by parame- 
ters which characterize a microstructure of polycrystals. 
In accordance with the above analysis, the thermal con- 
ductivity exhibits a crossover from a k ~ T 2 to n ~ T 3 at 
T ~ 0.1K for a chosen set of model parameters. It should 
be emphasized that such behavior of k(T) is specific to 
BWDD (which simulates a finite wall of edge disloca- 
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tions). For example, for the uniaxial dipoles one obtains 
that k ~ T~ 2 at low temperatures and k ~ for 
T — > 0£> , where 8d is the Debye temperature. 



In accordance with our scenario the mean free path is 
determined to be 
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FIG. 4. Reduced thermal conductivity due to biaxia 
WDD scattering, k x T -3 versus temperature T, calculatec 
according to (14) with the same parameter set as in Fig. 3 
Measured points for the, boundary-limited thermal conduc 
tivity in LiF (from Ref.E3) are indicated by triangles. 



III. BWDD IN DIELECTRIC GLASSES: 
THERMAL CONDUCTIVITY 

There are two important consequences of the previ 
ous section. First, it was found that the grain-boundary 
induced contribution to the thermal conductivity behave; 
like T 2 at very low temperatures. As is well knownE3, thi: 
behavior is peculiar to dielectric glasses, where k ~ T' 
for T < IK. Second, the critical temperature T de 
pends considerably on a size of the grain. This follow: 
from the condition 2qL ~ 1 which can be rewritten a: 
T w hv s /2Lka—wheve ks is the Boltzman constant (sei 
details in Ref.E^I). It particular, one obtains T ~ IK a 
2L ~ 20A. Hence a T 2 dependence of k can be extendec 
up to IK for materials consisting of microclusters witl 
average diameters of order 15 ~ 30A. It is intriguing tha 
exactly the same values are expected-,within the cluste: 
model proposed for dielectric glassesc3. 

This finding have stimulated— a detailed study of th< 
problem. It was shown in RcfJ13 that the experimenta 
data for the thermal conductivity in vitreous silica (a 
SiC^) can be explained by a combination of two scatter- 
ing processes. The first one comes from the sound waves 
scattering due to BWDD while the second one is the 
known Rayleigh type scattering which appear due to the 
local variations in structure. In this section, we present 
the model and extend an analysis to other glasses. 



m = (i- i 1 + i- t 1 ruc y\ 



(17) 



where lu comes from Eq.(15), and l s truc should be taken 
in the most general form describing the Rayleigh scatter- 
ing over the complete frequency range. The interpolation 
formula reads (see, e.g., Ref.o) 



lstruc — Y 1 (- ) 4 + Iq, 



(18) 



where Y is a constant which has been considered as a 
fitting parameter, and appears as the high-frequency 
limit. Fig. 5 shows lu, l s truc and 1{uj,T) with the model 
parameters for a-SiC>2 taken from Table I. 
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FIG. 5. Phonon mean free paths la (dashed line), htruc 
(dotted line) and l(u>) (solid line) as functions of frequency. 
The parameter set for a-Si02 is shown in Table I, K=2.6. 

One can see that l{ui) has a form typical for glassy ma- 
terials. At low frequencies, u) < 10 12 sec -1 , l(ui) ~ cj -1 , 
and the main contribution is due to the BWDD-induced 
scattering. In the intermediate region both scatter- 
ing processes are involved while at high frequencies the 
Rayleigh scattering becomes dominant. Notice that the 
region 10 12 sec _1 < w < 10 13 sec -1 is responsible for the 
plateau in the thermal conductivity. We have foundE£l 
that the size of this region decreases with 2L and/or Iq 
increasing. 

It is interesting to note that Eq.(17) supports the em- 
pirical relation l/X ~ 150, with A being the wavelength 
of a phonon, which holds for many glasses at low temper- 
aturesEI Indeed, at low frequencies / ~ lu- Spreading 
out Eq.(15) at qL << 1 one gets 



1 



A 2irD 2 (2vL)h 



(19) 



G 



This is a constant which depends on the model parame- 
ters which characterize the structural and elastic proper- 
ties of a material. It is reasonable to assume that these 
parameters vary only slightly in different amorphous di- 
electrics (see also Table I). 

This can explain the observed constant-like behavior 
of l/X. In particular, for our choice of parameters for 
a-SiC>2 one gets lu/X ~ 135. 

To calculate k with l(u>) from Eq.(17), it is convenient 
to use the dimensionless form of Eq.(14) 

1,4 T 3 r e D /T 
K = T^F^ / ~ l)- 2 l(x)dx, (20) 



where x = hu/kBT, and the specific heat capacity is cho- 
sen in the standard Debye form. The results are shown 
in Fig. 6 and Fig. 7. 

As is seen, there is a good agreement with the experi- 
mental data over a wide temperature range. Notice that 
we did not use any special fitting programs to get the 
best fit. Instead, we have fixed the parameters related 
to BWDD: a dipole separation 2L = 20 A, the density of 
defects n, = 2 x 10 11 cm~ 2 and the Frank index v = 0.1 
(except those for PS, see Table I), and tried to bring 
the parameters for Rayleigh scattering close to those in 
Ref.EJ. In our opinion, this provides better insight into 
the essence of the proposed model. Let us stress once 
again that a characteristic length 20A corresponds tp-the 
expected average size of clusters suggested in glassesE2tc3. 
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FIG. 6. Thermal conductivity vs temperature for a-SiC>2 
calculated according to Eq.(20) with l{ui) from Eq.(17) with-a 
set of parameters from Table I. Experimental data from RefH3 
are indicated by triangles. 




T( K) 

FIG. 7. Thermal conductivity vs temperature for a-Ge02 
(squares); a-Se (triangles) and polysfaa-ene (PS) (circles) (ex- 
perimental data are taken from Ref.EH). Theoretically calcu- 
lated curves are represented by solid lines. The set of the 
model parameters is given in Table I. 

IV. DISCUSSION 

The results obtained in the previous section call for 
an additional discussion. First of all, the main ques- 
tion arises of whether the proposed model captures the 
essence of the glassy state or a good fitting of experiments 
is purely accidental. At the moment we cannot provide 
the ultimate answer to this question. It is interesting, 
however, to consider the peculiarities of our approach in 
relation to another known ways of tackling the problem. 

As is well known, the very successful explanation of the 
very-low-temperature (T < IK) behavior of amorphous 
dielectrics has been done within the phenomenological 
TLS modelE3'Ej. Namely, at low temperatures the prin- 
cipal scatterers of acoustic phonons in glasses were pro- 
posed to be the tunneling states. As is known, however, 
this view posed some important questions. . First, the mi- 
croscopic basis for the TLS is still unclearEJH. Second, 
the quantitative universality seen in various glasses at low 
temperatures (for example, the above-mentioned relation 
l/X ~ 150) has not been explained. Namely, according 
to TLS model 

l/X cc (P)- 1 ex T g /V F , (21) 

where P is the density of TLS, T g is the glass transi- 
tion temperature, and Vp is the free volume frozen into 
the glass. However, the experimental data show that the 
above relation does.not depend explicitly on T g (see, e.g., 
discussion in Ref.Ell). 

Third, the universal properties of glasses in the 
intermediate-temperature range, 1 < T < 10 K, 
were viat— imdcrstood even qualitatively within the TLS 
modelE3 , E2lE3. Some of these problems, however, have 
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been solved later by invoking additional to the TLS con- 
cepts. In particular, a reasonable expression for the total 
phonon mean free path that allows to describe. fz, over a 
wide temperature range reads (see, e.g., RefsHE3) 



l(u,T) 



(22) 



where It and In are due to TLS and the Rayleigh scatter- 
ing, respectively, l a dd comes from some additional scat- 
tering mechanisms. As a possible candidate there was 
considered the .phonon scattering from some kind of dis- 
order (clusters&il, fractalstl, etc.). The modern theo- 
ries involve the ph Qn o n .| Sc attering from localized low- 
frequency vibrationsEZrE2E3 manifested themsekps both 
experimentally and by computer simulationscjO. All 
these successful approaches are essentially based on the 
TLS picture. Unfortunately, the nature of either TLS or 
the localized vibrations remains still unclear. 

Let us discuss these points in the context of our model. 
Our consideration assumes the new principal scatterer in 
amorphous dielectrics which has a clear physical origin. 
Indeed, as indicated above, there is a direct analogy be- 
tween BWDDs, finite dislocation walls, and grain bound- 
aries. Thus, the BWDD-based model-supports the clus- 
ter picture of amorphous dielectricsE2fL2 where bound- 
aries between clusters give rise to defects such as BWDD. 
In this connection the question arises: is there any sound 
experimental evidence for the existence of either clusters 
or BWDDs in glasses? At present the answer is no. As a 
possible reason one can mention a too small average size 
of clusters as well as their random distribution. In such 
an event, both clusters and BWDD are difficult to detect. 
Let us stress that we have considered in our model the 
randomly oriented dipoles (see averaging in Eq.(4)). 

In metals, however, the phase with long-range orienta- 
tional order and no translational symmetry ha 
perimentally observed by the x-ray scattering)! 
finding confirms the proposed point of view that super- 
cooled liquids and metallic glasses can be viewed as de- 
fected states, includi ng , cl isclinations, of icosahedral bond 
orientational orderBrir 4 !. In particular, in two dimen- 
sions liquids are regarded as hexatic fluids interrupted 
by point disclinations (i.e. local points of 5- and r 
symmetry)LJ. In accordance with their scenari 
there is a two-stage pairing process: disclinations first 
pair to form 5-7 dipoles regarded as dislocations which 
then pair at lower temperature to form a crystalline solid. 
It is interesting that grain boundaries are suggested to be 
linear arrangements of the form -5-7-5-7-5-7-. 

Regarding Eq.(21), we have shown in the previous sec- 
tion that the WDD-based model predicts a constant for 
the relation l/X at low temperatures (see Eq.(19)). It is 
interesting that this constant depends only on the model 
parameters which characterize the structural properties 
of amorphous dielectrics. Let us stress also that though 
Eq.(17) takes into account only two principal scatterers, 
it allows to describe the thermal conductivity of various 
dielectric glasses over a wide temperature range. 



This 



Another important question is the Law-temperature 
specific heat of glasses, C v . As is knownEH, C v is charac- 
terized by an anomalous linear temperature behavior. An 
explanation based on the TLS model looks quite correct. 
While this problem is beyond the scope of our paper, nev- 
ertheless, we can discuss briefly an expected contribution 
to the specific heat due to BW-D.D.. A similar problem has 
been considered long time agcO'tj. Granato analyzed the 
pinned-dislocation contribution^^) the specific heat and 
found that at low temperaturescH 



C v — 



pn' 



rida 2 Nks 



z e, 



T, 



(23) 



where p — v s *J p/G (G is the shear modulus and p the 
density of a material), rid is the dislocation density, a 
the lattice constant, Z the number of atoms per unit 
cell, and N the number of atoms per mole. This result 
has been, discussed in connection with dielectric glasses 
in Ref.El In particular, it was shown that there is a 
satisfactory agreement with the experimentally observed 
data for some glasses. As is knownt£l, with decreasing 
of the dipole separation 2L the biaxial WDD becomes 
equivalent to an edge dislocation with the Burgers vector 
b = 4Ltan(Q/2). Thus, the above result in Eq.(23) for 
dislocations should be valid for the BWDD-based model 
where 2L is expected to be very small. 

Notice also that a possible explanation of the specific 
heat behavior both at T < IK and 1 < T < 1£LK has been 
given within the soft atomic potential modeled as well as 
within the elastic dipole modelE3. These approaches in- 
terpret the specific heat peculiarities in terms of the TLS 
states and the additional quasilocal harmonic modes. As 
we have already mentioned, the excess harmonic modes 
coexisting with sound w aves .below ITHz actually have 
been observed in glassesc2lcll. It should be recognized 
that the presence of localized harmonic modes is pecu- 
liar to elastic materials with extended defects as well. 
In particular, the phonon spectrum in the presen ce ,, o f 
a dislocation was shown to possess localized modesElEl. 
Besides, the effect of localized vibration modes due to lin- 
ear defects on the thermal properties was studied within 
the fr amew ork of the vibrating string model of a dislo- 
cationEDE3. It is clear, however, that this consideration 
should be accompanied by a detailed study within the 
BWDD model. This problem invites further investiga- 
tion. 

Finally, let us discuss briefly a possible universality 
of the BWDD-based picture in polycrystals and glasses. 
In accordance with our results, a principal feature that 
distinguishes polycrystals from glasses is the size of a 
cluster. For example, let us consider the alkali halides 
LiF, NaCl, and KBr:KCN. The first two are polycrys- 
tals with the corresponding transport characteristics (see 
Sec.nl while KBr:KCN is an example of the orientational 
glassEl Within our scenario this markedly different be- 
havior can be explained by the lower average size of a 
cluster in KBr:KCN. As we have shown above, for this 
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reason alone the calculated crossover temperature can 
be moved from very low temperatures to IK. This means 
that, e.g., at 2L ~ 1000A the Independence of n will 
appear above O.IK. In addition, one can expect that for 
such large clusters the Rayleigh scattering will be sup- 
pressed while the umklapp processes become of impor- 
tance at high temperatures. As a result, the behavior 
of k will be typical for polycrystals. On the other hand, 
for small clusters the T 2 dependence of k goes up to 
IK, thenj-the Rayleigh scattering becomes important. As 
is knownd, there are no umklapp processes for randomly 
distributed small clusters. Thus, the proposed within our 
model scatterers (BWDD and the Rayleigh-type) will be 
of importance even at high temperatures. It is interest- 
ing to note that while the Rayleigh scattering becomes 
dominant at high temperatures, the proper behavior of k 
can be obtained only providing that the BWDD-induced 
phonon mean free path tends to a constant. 



ACKNOWLEDGMENTS 

This work has been supported by the Russian Founda- 
tion for Basic Research under grant No. 97-02-16623. 

APPENDIX A 

Let us find the exact expression for the perturbation en- 
ergy in Eq.(l). The WDD-induced strains Eij can be 
found by using of the Hooke's law 



(1 + a)4 - ao$5i. 



(Al) 



where afj are the stresses due to the WDD; /z and a are 
the shear modulus and the Poisson constant, respectively. 

For the chosen in Sec.LLgcomctry (sec Fig.l) the WDD- 
induced stresses of,- areO 

l£l [1 (x + L) 2 + y 2 | y 2 



ri. 

2tt(1 -a) 12' 



(x - L) 2 + y 



(x-L) 2 + y 2 (x + L) 2 +y 2 



-h 



{x + L)({x + L) 2 -y 2 )) 



+ k 



(x + V) 2 + y 2 
{x-L)([x-L) 2 -y 2 )) 



{{x-L) 2 + y 2 ) 



(A2) 



+ 1. 



(x-L)((x-L) 2 + 3y 2 )) 



((x-L) 



2 _|_ y2 



(A3) 



x + L 



tt(1-(7)L2 (x-L) 2 +y 2 L (x + L) 2 + y 2 

x — L 



+ k 



(x - L) 2 + y 2 



(A4) 



Notice that Eqs.(A2)-(A4) describe the stresses due to a 
finite wall of edge dislocations at large distances. Apply- 
ing Eqs.(A2)-(A4) in Eq.(Al) one gets all the components 
of the strain tensor Eij and, finally, Eq.(2). 

APPENDIX B 

The perturbation energy given by Eq.(2) takes the fol- 
lowing form in polar coordinates 



U(r, 4>) = B 



1 r 2 + 2rL cos (f> + L 2 

2 r 2 - 2rL cos 6 + L 2 



h- 



r cos 6 + L 



+ 1 



2 + 2rL cos 6 + L 2 r 2 - 2rL cos 6 + L 2 



(Bl) 

The matrix element in Eq.(5) with the perturbation en- 
ergy from Eq.(Bl) can be calculated using the following 
formulas: 



00 „2fc-l 



E 

k=l 



cos(2fc - 1) 



2k- 1 



1 1 + 2zcos i 
- In 

4 1 — 2zcos ( 



z z < 1, 



(B2) 



z k cos ktf) = — 



1 — Z COS ( 



fc=0 



— 2z cos 6 + z- 



r, |«|<1 (B3) 



Substituting Eqs.(B2) and (B3) into Eq.(Bl) and inte- 
grating in Eq.(5) one obtains 



U (p, a) 



AttBL 



pS 



Jo{pL) cos a 



ufi rl, (x + L) 2 +y 2 (x + L) 2 
— in ■ 



27r(l-cr)L2 {x-L) 2 +y 2 (x + L) 2 + y 2 



+ 



^(-l) fc J 2fe (pL)( 



cos(2fc + l)a cos(2fc — l)a 



fc=i 



2k + 1 



2k - 1 



{x _ L) 2 (x + L)((x + L) 2 + 3y 2 )) 

— h — 



(x-L) 2 + y 2 



{x + L) 2 + y 2 



2ttB a 
pb 



JaipL) cosa + 2cosa^(— 1)* J2k(pL) cos2fcc 

(B4) 



k=l 
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where A; = li — 12 and J m (z) is the Bessel function. The 
first term in Eq. (B4) comes from the integration of loga- takes the form 
rithmic function in Eq.(Bl) while the second one comes 
from two last terms in the r.h.s. of Eq.(Bl). ^-1 _ 

We have used the following standard integrals deriving h 2 v 2 | 2q 3 

Eq.(B4) 



The exact expression for the mean free path in Eq.(3) 



n,B 2 ir 1 



A, 



— h(qL) + ^h{qL) - -ih(qL) , 



2q 



(Bll) 



2tt 



J exp(iz cos 4>) cosm<j)d(j) — 2Tri m J m (z) , (B5) 



where 

h(qL) 



dO 



-(1 - J (4 9 Lsin(6>/2)) 



J z k J k ^ 1 (z)dz = z k J k (z). 



(B6) 



The first sum in Eq.(B4) can be simplified by differen- 
tiation with respect to a. The second sum in Eq.(B4) 
reads 

00 11 
^(-l) fc cos(fca)J 2fe (z) = -cos[zcos(a/2)] - -J (z) 

fe=i 

(B7) 

After straightforward calculations one obtains 
B ■ 

si~ p 



sin 2 (0/2)' 

l6q 2 L 2 Tr(j 2 (2qL) + Jl(2qL)j - SqLirJo^qtyJ^qL), 

(B12) 



/■ 2 vi 

h(qL) = J (-+Jo(4gLsin(0/2))- 



Ji(4gZ,sin(0/2)) 



AqL sm{0/2) 

^(l + Jo(2«i)) - ir(j 2 (2qL) + J 2 (W), (B13) 



dO 



U(p, a) 



Am 2-KiAi 
— ^- sm(pL cos a) H cos a cos{pL cos a) 



P 



r 2TT 

ML) = \ 
Jo 



de 



Ji(4 9 Lsin(6>/2)) 



To find 5ft(0) in Eq.(4) one has to make averaging of 
\U(p, a) 1 2 over a: 

1 f 2n 

\U(p)\ 2 = |< q|C/(r)|q' >| 2 = —J \U(p,a)\ 2 da = 



sin(6>/2) 
4qLir(j 2 (2qL) + J 2 (2qL)y 



(B14) 



The final result for the mean free path is given by 
Eq.(7). 



2irB 2 
S 2 



C(7 



2 A 

sin 2 (pL cos a) 5- cos a sin(2pL cos a) 



A 2 \ 

H 5- cos 2 a cos 2 (pL cos a) )da. (B8) 

p z ) 



Using 

n 



cos(zcosa)\ , . , „ /cos(ri7r/2)\ , 
sin(zcos a) ) C ° S{na)da = 2 %in{nn/2)) J ^ 

(B9) 



one finally gets 

4tt 2 B 2 f 2 . ,A 2A ; 



\U{ P )\ = 



S 2 p 4 



(l-J {2pL))--^J 1 {2pL) 



2pL 



Substituting Eq.(BlO) into Eq.(4) one obtains the effec- 
tive differential scattering radius in Eq.(6). 
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TABLE I. Parameters used for the numerical fits. Units 
are the following: v s (10 5 cm/sec), 6d (K), nt (10 11 cm -2 ), 
2L (10~ 7 cm), Y (cm-'K- 4 ), l (10~ 7 cm). 



Material 


V 3 


Qd 




2L 


V 


lo 


Y 


a-Si0 2 


4.1 


342 


2 


2 


0.1 


1.5 


1 


a-Ge0 2 


2.6 


192 


2 


2 


0.1 


0.6 


2.9 


a-Se 


1.19 


113 


2 


2 


0.1 


0.2 


90 


PS 


1.67 


123 


5 


2.4 


0.1 


0.5 


80 
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